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Abstract 

In this paper we deduce a quite general formula which allows the relation of clock rates at two different 
space time points to be discussed. In the case of a perturbed Robertson- Walker metric, our analysis leads 
to an equation for the comparison of clock rates at different cosmic space time points, which includes the 
Hubble redshift, the Doppler effect, the gravitational redshift and the Rees-Sciama effects. In the case of 
the solar system, when the 2PN metric is substituted into the general formula, the comparison of the clock 
rates on both the earth and a space station could be made. It might be useful for the discussion on the 
precise measurements on future ACES and ASTROD. 
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1 Introduction 



Recently, atomic clocks with a time-keeping accuracy of the order of 10~ 18 in fractional frequency have been 
considered ^ |2] • Also a spatial experiment named ACES (Atomic Clock Ensemble in Space) mission [1] is 
scheduled to be launched in near future by ESA (European Space Agency). The purpose of ACES is to obtain 
an accuracy of order 10~ 16 in fractional frequency. In such a situation (1CP 16 — 10~ 18 level), 2PN (second 
Post Newtonian) approximate framework has to be carried out before hand. Also, ASTROD (Astro-dynamical 
Space Test of Relativity using Optical Devices) 5 , 6 has been proposed which is now under consideration by 
the Chinese government. It is a three advanced drag-free spacecraft system, and mature laser interferometric 
ranging will be used. The accuracy of measuring 7 (about 10 -9 ) and other parameters will depend on the 
stability of the lasers or clocks. In ASTROD the desired accuracy of the optical clock is about 1CP 17 . This plan 
also needs a 2PN level on the comparison of clock rates and equations of motion for planets. The precision of 
2PN level on the comparison of clock rates (or time transfer) has been discussed in Ref. by means of world 
function. But as we know, the calculation of the world function is not easy [Hj • Also they discuss the frequency 
shift to the order of c -4 only in the field of an axisymmetric rotating body, not in a general case. Therefore 
we deduce a general formula in a different way. Our general formula can also be applied to cosmos and easily 
extended to an even higher order (higher than 2PN level). But here we have to emphasize that the purpose of 
our general formula is to discuss the frequency shift (or variation of the period) between two arbitrary space 
time points in a metric theory. It is not to attempt to establish the synchronization of the distant clocks, e.g. 
a grid of admissible radar 4-coordinates [UEII; which is very useful for a GPS (Global Positioning System) 
system of spacecrafts. It is limited only to the clock rate, not the synchronization of the clock. It does not 
involve the splitting of the space time. 

We hope our general formula might be useful in ACES, ASTROD and LISA (Laser Interferometer Space 
Antenna) [llj which is a worldwide interested space mission. Our formula might be also useful for the precise 
measurement of Hannay effect ^JElEl (the Doppler shift variation is equal to 3.7 x 10 -15 for a geostationary 
satellite). Hannay effect is the classical counterpart of the quantum geometric phase discovered by Berry |15|. 
In future experiment, we have to know about the higher precision of the modelling frequency shift before hand. 
In this way, our formula may be useful for the establishment of such a model. 

On the other hand, one of the most basic experiments in physics is the measurement of time, often at different 
positions in space. Many formulae have been suggested for the comparison of clock rates at different positions, 
based on certain simplifying assumptions about which effects are dominant. We mention the most prominent 
ones. The change of the clock rates can be related to: the relativistic Doppler effect, the gravitational redshift, 
the Hubble redshift, the Rees-Sciama effect and so on. Let us recall these simple formulae on the comparison 
of clock rates in the following. 
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In special relativity, the Doppler effect is shown by 

a 1 — VAB ■ -h- 

At b _ 1 c |fc| m 

where Ata and Atb are the proper period of a light signal emitted by a source A and received by B respectively, 
vab is the relative velocity from A to B, k is the wave vector in flat space from A to B. The relativistic Doppler 



effect (or transverse Doppler effect) 1/y (1 — ^y) was measured long time ago [T^]. Recently the time dilation 
factor has been measured with an accuracy of 2.2 x 1CP 7 |17| . 
The gravitational redshift is given by 



At b = y/-goo(B) ^ 
Ata ^J-g m {A) ' 

In Newtonian limitation the metric tensor goo = — 1 + 2w/c 2 , then Eq.J2J) becomes 

^ = 1-^ + ^ + 0(4), 
Ata c c 

where w is the Newtonian gravitational potential, 0(4) is the abbreviation symbol for 0(c -4 ) as well as 0(n) 
for 0(c~ n ) in the following formulae. The gravitational redshift was measured with an accuracy of 7 x 10 -5 

dHj. 

The Hubble redshift is written as 

Atb _ R{Bl . , 

Ata R{A) ' W 
where Atb and Ata are the period of receiving and emitting signal, R is the cosmic scalar factor. 

When photons pass freely through a region of voids or supercluster, an additional redshift (or blue-shift) will 

be created which has to be added to the normal Hubble redshift. Because of the pioneer work by Rees & Sciama 

|19|. this additional redshift due to the voids or superclusters is called as Rees-Sciama effect. In their original 

paper they modelled the voids or superclusters by a spherically symmetric density over a Swiss cheese embedded 

in the Robertson- Walker metric. The physical meaning of the effect is explained in detail, but no derivation 

from the first principle is given. Later many scientists |201 1211 12"2~] discussed the anisotropy of CMBR(Cosmic 

Microwave Background Radiation) caused by inhomogeneous Universe (Rees-Sciama effect). The original (Ress- 

Sciama) spherical Swiss cheese model has been extended to an integration of the time dependent Newtonian 

potential [221 

^Ata^ = ~ w{tBl ^ + w(tA ' X%a) + 2 J A ^ dt ' (4) 
where w is the Newtonian potential. 

The physical conditions causing all these effects may be all present. In early 90's the Hubble redshift, the 

gravitational redshift, the Doppler effect and the Ress-Sciama effect have been combined into one equation (see 

Eq.(6) in |24| ) in first order approximation: 

l + ^=|^y|l + ^(^-0 o ) + 2^ T °drl.V ( /» + n.(v e -v o )| , (5) 
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where z is redshift, <j> is gravitational potential, the subscript e(o) denote the emitting (observer) point , 1 = k/fe D 
(k a is the tangent vector to the null geodesic connecting the emitting point and observer), v is 3- velocity. Since 
all of terms are the level of the first order approximation, the coupling terms do not exist. Martinez- Gonzalez 
et al. [23] do not deduce the equation (JjjJ through an exact method, it is difficult for us to extend the formula 
to higher order precision. At last, the coefficient 5/3 may not be clear. 

Accordingly, a comprehensive approach, starting from the first principle, is needed in which the physical 
conditions for all these effects are taken into account at the same time. Such an approach should lead us 
to a synthetic formula which reflects all these effects in a compact way and which should provide additional 
information, due to possible interactions which could not be incorporated in the isolated approaches for the 
individual effects. 

Normally the clock rates of A and B are represented by the proper time interval Ata and At b at position 
A and B respectively. In general a comparison of the clock rates between Ata and Atb by means of differential 
coordinate time At a and At B in global coordinates can be achieved. The relation of the coordinate time between 
A and B is established by null geodesic line (light-ray). The differential coordinate time on null geodesic line 
can be easily calculated out 25 



The minus and plus sign are taken in I and III quadrants (in x-t coordinates) and II and IV quadrants respec- 
tively. Normally we take the minus sign. Using these assumptions we can obtain a general formula for the 
comparison of clock rates by means of some integral. 

In section II, we first time deduce a general formula by means of "calculus of differences". In section III, 
substituting the simplest perturbed Robertson- Walker metric into the general formula, we obtain a formula for 
the comparison of clock rates at different cosmic space time points, which includes the Hubble redshift, the 
Doppler effect, the gravitational redshift and the Rees-Sciama effects. In section IV, by using the 2PN metric in 
multiple coordinates the 2PN comparison of clock rates on both the earth and a space station in the solar 
system is made, being possibly useful for the precise measurement of ACES, ASTROD and LISA in future. 
Some conclusion are offered in section V. 

2 General formula 

In a global coordinates (ct,x l ), a source A moves with a velocity v A and a receiver B with a velocity v l B . The 
clock rates in A and B are directly related with their own proper time Ata and At b . To compare them, we 
need to know the relation between the coordinate time interval At a and Ats, because 



cdt = 



g 0i dx l ± ^(goigoj - go gtj)dx l dx3 



(6) 



5oo 



At a _ At a At A At B 
At b ~ At A At B At b ' 



(7) 
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Since — c 2 dr 2 = ds 2 , therefore if the velocity of a standard clock (A or B) in the global coordinates is v, we 
have 

At B = ATB (8) 

sj -[900(B) + 2g 0i (B)v B /c + gij (B)v B v B /(?} 

At A = ATA =, (9) 

V -[5oo(A) + 2goi{A)v\/c + g l3 {A)v\v\lc 2 \ 

where g^u (goo, 9oi and gij) are the global metric. As abbreviation we introduce 

Ga = -(.900(A) + 2g 0t (A)v l A /c + 9l] {A)v\v A /c 2 ) , 
G B = -(900(B) + 2g oi (B)v B /c + 9ij (B)v B v B /(?) . 

Here At, being an integrable coordinate time interval, has unique meaning throughout space. Eq.J7|) is known 
for long time. In most discussions on the gravitational spectral shift, people tend to consider a metric as 
independent of time (stationary metric) and think points A and B are fixed in space, then take At A — At B 
for granted |27j . One of the main purpose of our paper is to calculate the relation between At a and Ats by 
means of the "calculus of differences". Assuming that, at t Al (coordinate time) a source A emits a first pulse at 
position Ai(x Ai ), then a receiver B receives the first pulse at position B\(x l B ) at time ts ± - A second pulse is 
emitted from A at position Ai[x l A ^) at t^ 2 , which is received by B at position B-2.(x l B ^) at time tB 2 - Two pulses 
transmit along the null geodetic line 1 and 2 respectively. The emission time and reception time are related as 
follows: 



. . , 1 f Bl -90idx l - Vigoifjoj - goo9i )dx l dx3 

tj}\ = Mi + " ! (.J-U) 



c J At 9oo 



. 1 [ B2 -9wdx l - ^(goigoj - goo9i )dx l dx3 

tB 2 = M 2 + - / ■ (11) 



c Ja 2 9oo 

We define dx 2 = 5ijdx l dx\ the geometric meaning of dx is the spatial differential length of the line in the flat 
space. Then the relation between the emission time and reception time can be rewritten as 



I / " ~9oi^t ~ \J (9oi9oj ~ 9oo9ij)^^ 



dx 1 dxi 



t B =t A + - dx. (12) 

c J a 9oo 

In a weak field, goigoj is a small quantity (~ 0(6)) and the spatial conformal isotropic condition |281 [26 is 

9oo9ij = -Sa - ^ + °( 6 ) . ( 13 ) 
where q^j is a spatial anisotropic contribution in the second order. Then Ea. (|12fl simplifies to 



1 f B 1 



tB = t A / 

c J A 900 



dx 1 qtj dx 1 dx° 
dx 2 c 4 dx dx 



dx + Q(6). (14) 



Briefly, we define 



~9oi%- - V (9oi9oj - 9oo9ij)-^--^ 

F(t, x l ) EE 1 , 15 

C300 
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Eci. (|12[) becomes ts = + J A F(t,x l )dx. Considering the time interval of two pulses as a period of the light 
frequency (an atomic clock), it is very small (about 10 -15 sec), during which the source or the receiver move 
a distance only about a magnitude of angstrom (A). Therefore in the following treatment only the first order 
approximation (linear approximation) is taken (all of higher order terms are neglected). From now on, we take 
Ai as A, A 2 as A + AA and as well as B\ as B, B 2 as B + AB, At a — tA 2 — tA 1 , AtB = ts 2 — is^ According 
to the "calculus of differences" , we have 



where 







At B -- 








A 


/ F{t,x l )dx 






J A 





At, 



fx(B+AB) 



F{t,x l )dx 



F(t,x l ) 



dx 



lino 2 



■(B) 



F{t,x i 



x{A) 



dx . 



(16) 



(17) 



lincl 



c(A+AA) 

If the metric is independent of time and points A and B are fixed in space, Eci. (|16|) will lead to At a — Ate just 
as we mentioned earlier. Note that in Ea. l|16|) we are dealing with finite differences. 

The null geodetic line 1 and line 2 are very close, but different. Therefore we can rewrite the integral along 
the line 2 as 



x(B+AB) 



x(A+AA) 



F(t,x l ) 



dx = 



linc2 




F(t + At, x + Ax)dx , 



(18) 



rx(B+AB) 
L(B) 

where At and Ax" 1 means the deviation from line 1 to line 2 (certainly At and Ax 1 are the function of space- 
time point). Substituting Ea. (|I8(l in to Eq. l|17(l and considering a linear approximation in the "calculus of 
differences", the difference of the integral in Ea. (|16[l can be divided into three parts: 



Fdx 



x(B+AB) 



Fdx 



x{B) 



x(A+AA) 



Fdx 



x(B) 



AFdx, 



x(A) 



where 



x(B) 



AFdx = 



x(A) 



x(B) 



x{A) 



(F(t + At, x l + Ax 1 ) - F(t, x 1 )) dx . 



Ax(A) and Ax(B) are given by 

Ax (A) = x(A + AA) - x(A) 
Ax(B) = x(B + AB) - x{B) 



\k A \ ' 


G?X 

dx 


Ax = 

A 


\k A \ 


k B 


dx 


Ax = 

B 


k B 


\k B \ 


dx 


\ks\ 



VAAt A , 
■ v B AiB 



(19) 
(20) 

(21) 
(22) 



Here \<la is the wave vector at point A of the light signal emitted from A and k^ is the value of the wave 

^xby^ 



vector at the received point B. In Eas. (|21() and (1221 we have replaced ^§-Ax by ^§-Ai. Because At a is a small 



quantity in the right hand side of Ea. H21l) . there is no difference between k(Ai) and k(^2) at the level of linear 
approximation. We simply take k^ to substitute k(^4i) and k(^2), as well as in Eg. 1)2211 ks = k(£?i) = k(i?2). 
The first and second terms in Eg. l|19|) can be written as 

r-x{B)+x(AB) 
Jx(B) 



Fdx 



(23) 



G 



x(A)+x(AA) . 

Fdx = -F(A)^-^Ai A . (24) 

e(A) 

Eg. (|20fl is the difference of the integral between line 2 and line 1 when = vb — (i.e. the boundary of 
integral is fixed), which can be expanded as 

r x (B) r*{B) f at? f)p \ 

/ (F(t + At,x l + Ax l )-F(t,x l ))dx = / — A< + — Ax 1 ) dx . 



(25) 



As we know, if F is independent of time, then f*ffi AFdx = f*f$ j^iAx l dx = 0, since for fixed boundaries 
the light ray is unique (no deviation). When F is dependent on time, there are two curves. The second term 
Jx(A) §x T ^ xl< ^ x cause d by time dependent metric is a higher-order term compared with Q^-Atdx, i.e. 

X ( B ) dp r x ( B ) gp 



— Ax l dx< / —Atdx. (26) 

x(A) dx 1 J x[A) dt 

In Ref. [2H! it is also mentioned that ^-Aafdx does not cause the change of frequency, even for cosmic 

distance. Finally we substitute Eq.J2SJ, lEH, 123) i nto Ea.lfIS)l. and use Eq.ljSJ and Eq.© to represent At a 
and Ats by At a and Atb- The general formula is obtained as 

At b [Gb~ ( 1 - F ( A ) ) ^\ VGb~ r (B) fdF^,dF 



-, , . , , I , . — At+— Ax 4 dx. (27) 

Ar A V G A \l- F(B)^ J Ata i _ F(B)V^) V) V & dx> > 



\k B \ 

Eq. (|27|l is our main result. Based on this equation we can calculate the comparison of the clock rate between 
two arbitrary points of space time. The equation is dependent on the metric and the path of the null geodesic 
line. The metric can be solved from the field equation (the general relativity certainly is a most desired (but 
not the unique) theory) in different situations (the energy momentum distribution and boundary conditions) 
while the path of the null geodesic line can be obtained from null geodesic equation. We do not care about 
the individual solution, even in the following examples we only want to show how the formula is efficient in 
deducing some rather complex formula (such as Rees-Sciama effect) smoothly (see Section 3) or the coupling 
terms (see Section 4). 

As an example, we consider the Doppler effect of a moving source in Minkowski metric 



-1 
6 



since v B = 0, G B = 1, G A = 1 - and ^ = Jg- = 0, F(A) = ±, so that 



1-^ 

This is just the formula of the Doppler effect in the special relativity (Eq.QJ) 



At B c \k A \ , , 

~^--i=k- (28) 
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The other simple example is the gravitational redshift. Considering a static gravitational field (e.g. Schwarzschild 
metric), in which both source and receiver without moving (v^ = vg =0), the general form then becomes 



At a y -.goo (A) c 2 c? 

where the last step of above equation is the Newtonian limitation. Ea. l|29|l is just the equation shown in normal 
textbooks of gravity (see Eq.©). 

Naturally, for an explicit evaluation of Equation (|27f) . the light-ray trajectory should be known. If we take 
dl = cdt as the arc length of light-ray trajectory, then we get [3"U| 

dx l , dx\ v dx\ 

-dl= n+a ^ +a ^ + --> <™ 
where n l is a constant unit vector, parameters ot\ and ai ■ ■ ■ are equal to zero or 1 dependent on which 
approximate level we consider. For 1PN level [a\ = 1, a» = when % ^ 1), we have |31| 

n-^ = -^, (31) 
at c 



= 2\Vw- 2n(n • Vw)] , (32) 



dt 2 

For 2PN level (cci = a% = 1, a, = when i > 3), by means of iterative method we deduce 




f|nW , (33) 



dx\ p \ ( 2w tj n j w t 2w (j,k) i k 2w.jdx[ p \ 

! ~~ ~ ^2—" 71 + c 2 dt J 

c 2<r 

where u>, io' and qij are defined in |281 I26| . Hereafter we use the abbreviation DSX for "the papers written by 
Damour, Soffel and Xu" H3 E2 . In principle, if we knew the metric to N-PN level, we could solve 
(or up to such a N approximate level. Therefore =%- can be taken as a known function. 

In cosmological problems, dx l has to be replaced by R(t)dx l where R(t) is the cosmic scalar factor (or the 
expanding factor) |35|. x l is a comoving coordinate now. Then Ea. (|30|) has to be written as 

dx 1 dx]„ dx%„ 

R(t)=r = n l + a x R{t)—^- + a 2 R(t)-^ + ■■■ . (35) 
dl dl dl 

Usually, in most problems of cosmology, the deviation from the leading term would not be considered. The 

relation between dx and dl can be simply obtained from Ea. l|3U[) for ot\ = oti = 1 

UxV = i + 2n dj^ + 2n dj^ + dj^ d^ L + 
\dl ) dl dl dl dl 

Now we shall discuss two applied examples in the following section. 
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3 Application in Cosmos with perturbed R-W Metric 



First we recall the unperturbed Robertson- Walker metric 



dg2 = _ cHe m%jdM 

(l + fr 2 ) 2 



where R(t) is the cosmic scalar factor, k = — 1, 0, +1 is corresponding to the open, flat and closed universe 
respectively. R(t) has the dimension of length and dx % is dimensionless. As we already know R(t) is model 
dependent. In cosmology, by using Ea. l|37l) the light ray ds 2 = yields 

dt \J 5ijdx l dxi 



R{t) 1 + \r 2 
The right side of Eq.(j2EJl is independent of time, therefore 

dt dtA dt b 



(38) 



(39) 



R{t) R{t A ) R(t B ) ■ 

Since usually we do not consider the local gravitational redshift and the Doppler effect in the problem of 
cosmological expansion, we then have At a — Ata and Ats = Atb, and thus the formula for the Hubble 
redshift is 

Ata_ _ A^4 _ R{t A ) 

At b At B R(t B )- ( > 

The results of Ea. (|40|l can also be deduced directly from the general formula Ea. (|27|) . if we take va = vb = 

WA — WB = 0. 

Next, we consider a linearly simplest perturbed Robertson- Walker metric of the form 

,9 9 /„ 2uA , 9 / 2w\ R 2 5 il dx l dx j ,, , 

ds 2 = -c 2 1 - — )dt 2 + 1 + — !2 T , 41 

V c 2 J V c 2 J (l + | r 2) 2 

where the gravitational potential w = w(t,x l ) is assumed to be a small quantity. Certainly, we might choose 
more general perturbation, e.g. Sachs- Wolfe integrated effect or the gauge invariant perturbation |37U38| . In 
fact if you know any kind of metric, just substituting it into the general formula (Eq.J5JJ) to get the comparison 
of clock rates between any two arbitrary points. If we use some complex perturbation, the formula on the 
comparison of clock rates would be too long to see our extension clearly. By the way, the application of our 
general formula in discussing CMBR |39j might be very attractive for us in the future. For the moment we 
will concentrate on solving the comparison of clock rates between two arbitrary points. Here we consider the 
Doppler effect caused only by the motion of the source, then vb = (also possible va — 0, then vb ^ 0). The 
velocity of the source A is 

v A = R(t A )^. (42) 
F(A), Ga, and G B can be calculated as follows 



= [ - 7 V , kr l\ dX dX =^- + 0(3) , (43) 
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where we have neglected all of higher order terms and consider bijn % v? = 1 and kr 2 as higher order term also. 



The quantities G a and G B read as 

\/~Gb = \/-goo(B) = 1 - \w(t B ,x l B ) 



G , 



IG 



"<7oo — %9oi — — 9ij 



dx 1 



dx 1 dx 3 
~dt~dt 



1 



-w{t A ,x A ) 



~ V A 



1 



1 



TT = 1 + ~oW(t A , x\) - —w(t B ,x B ) + -f 



(44) 
(45) 
(46) 



Ga '" C ' ' 

Now we calculate the last term of Ea. (l27[l . We should pay attention to that, in Eas.l|44[l and 145(1 we use the 
same symbol w to represent the gravitational potential at point A and point B, but really their form maybe 
very different (along the right-ray a potential may caused by supercluster or the sun) , normally we can not use 
the same function to describe the potential in different points A and B. In the last step we consider the integral 
in Eq.JHJ) 



—Atdx , 



x(A) 



dt 



(47) 



where we have omit the term of f^ffi §§rAx' l dx because of Ea. i|2t)[l . Therefore we have 



dt 



(A) 



R(t) 1 



2w 
~7s 



AtW 6ijdx l dx 3 
c(l + fr 2 ) 



(48) 



In fact we know the relation between At a and At B , therefore at any time t, in first order of approximation 
from Ea. (|39|l we have 



Ea.QBfy becomes 



X{B) R{t)At A 

(A) R{tA) 



• /n 2uj , d(2w/c 2 )\ J8 l1 dx t dx3 
R(t) + R(t) — + R(t)- K ' ' » V J 



dt 



(1 



V) 



Considering null geodetic line, Ea. ((41|l yields 



(-5) 



dt = ± 



^ w \ R(t)y / Sijdx l dx : > 



C 2J 



(1 + !' 



which we use to evaluate the integral / and get 
At A 



I = 



R(t A ) J x{A ) 



dt = At a 



c 2 dt 

The second term of Eq. (|27|) then becomes 

fGj f R(t B ) - R(t A ) 
R(t A ) 



R{t B ) - R(t A ) 
R{t A ) 



c 2 R(t A ) J x(A) 



R(t)-dt 



(49) 



(50) 



(51) 



(52) 



c 2 R(t A ) J x{A) 



r {B) n ^9w^ 

Jx(A) Ct 1 



(53) 
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Substituting Eas. l43Jl . l|46|l and (|53JI into Ea. l27|) (pay attention to that, in cosmology the term ^ • v should 



be replaced by & ■ in Eq.|(2TJ), (JUJ), and Ea.l|77|l1. we finally obtain 



At b 
At a 



(w(t A ,x A ) - w(t B ,x l B )) 



R(t B ) k A • v A 



(54) 



^R(t A ) c\k A \ c 2 R{t A ) J X ( A) 

where ^ (w{t A ,x A ) — w{t B ^x l B )^ is the contribution from the normal gravitational redshift; ~nffi and ^ are 
the Doppler effect and transverse Doppler effect (or relativistic Doppler effect) respectively; R(t B ) / R(t A ) just 
contributes to Hubble redshift; the last term is related to Rees-Sciama effect. Ea . I|54[) can be compared with 
Eq.Q if we put va = and R(t) = R(t A ) = R(t B ), then 



1 2 f xlyB) dw 

1 + -x (w(t A , x\) - w{t b , x B )) + 3 / — dt , 

c c Jx(A) 01 



(55) 



where w is the same as u in Ref. |23j where c = 1 units is taken, and their results totally agree with ours. 
We thought that our general formula allows to derive the Birkinshaw-Gull effect 0H1 QD t° ' if we use suitable 
perturbation functions for w(t,x z ) and Wj(t,x l ). This will be discussed in another paper. 



4 Application in Solar System with DSX metric 

Our formula l|27|) is of substantial generality and accordingly allows a great number of fruitful applications. 
Some of these have been given in the previous sections. However a important application probably is to the 
solar system since there we have a chance to do measurement and thus by comparison with measurements our 
formula can be tested. In this section, we provide a fundamental formula with the precision 0(4), but we have 
not expanded the potential and vector potential by means of multiple moments that has to be done in the 
practical problems. 

In near future high precision measurement will be done up to 2PN level as we mentioned before, thus 
allowing coupling terms (i.e., the terms connecting the gravitational redshift, the Doppler redshift and so on) 
to be measured. Our scheme (the general form Eq. (1271) ) offers the possibility for this if an appropriate 
assumptions about the metric are used. Accordingly, in this section we start from DSX formalism |281 1321 [33] 
and its extension and evaluate formula (|2*7jl for this metric. This extended DSX metric is described by 

.g o = -cxp^-^^ +0(6), (56) 
Am- 

g^ = — r + °( 5 ) . ( 5? ) 

. 9y -=%exp(^f) + |f + 0(6), (58) 

9ij9m = Sij ~ |f + °( 6 ) ■ ( 59 ) 

In fact, in the following calculation, q^j appear only in the function F(x l ,t), but in the final 2PN formula of 
clock rates (to see Eq.{2J0 qij does not exist which agree with the result in Ref. 0. Substituting Ea. l(5^ - (|5§|) 
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into Ea. l27|) . calculating all of components, we can get a general formula for the comparison of clock rates at 
2PN level. 

We begin by evaluating the terms Ga and Gb for the extended DSX metric: 



Ga = -900(A) -2g 0i (A)^- gij (A)^ 



. 2w(A) + 2w 2 {A) + 8wi(A)v A v\ 2w{A)v\ 



2w{B) 2w 2 (B) 8 Wl (By B 2w(B)v 2 B 

^ B -L n > A * A O A 



0(6), 

1-0(6). 



(60) 
(61) 



Next we consider the terms F(A)\s.a • v^/lfc^l and F(B)'k.B ■ vs/|fcs|. Since • / 1 A; a | and k^ • v.b/|&b| are 
first order already, so F(A) and F(B) need to be calculated up to c -5 level. From Ea. i|15fl we get 



F(A) = 



CffOO 



dx' 




dxi 




J-googij ^ 


dx 



cgoo 



where we have neglected ffoiSojC^ 0(6)). By using Eq. 



Ea. (|62|) then becomes 



F{A) 



and similarly 



F{B) 



2w(A) 4wi(A) dx 1 
c 2 c 3 dx 



2w(B) 4wi(B) dx* 
c 2 c 3 dx 



2w 2 (A) 1 
c 4 + 2? 



2w 2 (B) 1 

C 4 + 2^4 



dx dx 



dx 1 dx 3 
dx dx 



0(6). 



0(6) 



(62) 



(63) 



(64) 



At last, we consider the integral (the second term) in Ea. (|27|l 



B 'dF A dF . , 

At + — Ax 1 dx. 



dt~ 1 a,-- ' (65) 

where because of Ea. l|26(l we will omit the second term in above equation. Accordingly instead of (|65() we 



consider 



dF 



Atdx. 



(66) 



1 A dt 

As we know At a at A and Ats at B, we could calculate At at an arbitrary point between A and B approximately. 
Then the definite integral l|66l) can be evaluated by the median method, i.e., 



^ dF , / >s dF 

—Atdx = At / — dx . 



dt 



dt 



where At is the median value, for which we introduce a parameter 77: 

At = r]At A ■ 



r\ is a value closed to 1. For the term ^ one finds 



dF 1/2 dw 4 dwi dx 1 Aw dw 1 dqij dx 1 dx 3 



dt c\c 2 dt c 3 dt dx c 4 dt 2c 4 dt dx dx 



0(6). 



(67) 



(68) 



(69) 
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In the solar system the change of the potential (the other metric much smaller) with time is very small (less 
then 0{2) level), only the leading term is considered. Therefore the term ^ ^fdt is already on 0(4) level, 
but not 0{2). Furthermore, if we take Ea. (|3l)|) into account, we have dx = dl + 0(2) = cdt + 0(2). Then Ea. l(6"9"|) 
simplifies to 



OF , 2 
-ax = -7T 



dw 
~dt 



dt + 0(6) . 



i a dt c 

If we consider quick variable field (e.g. field in pulsar) we have to take Ea. i|t)9[l to substitute into Ea. (|70|) . 
The second integral term of Eq. I|27|l then becomes 

r-B 



(70) 



2 V 



<(Ja 



dw , 
-rrrdt. 
k s-VB \ „2 J a at 



(71) 



Gathering all evaluations done thus far in this section we arrive at following general formula for the solar system 



Ars 

AT4 




1 - F(A) 



k-A ■ va 2r? 

\k A \ c2 



dw 
~dt 



dt 



(72) 



In Eq. l72|l we only consider the leading term, namely the scalar potential changing with the time, but in our 
scheme all of the second post Newtonian terms (see Ea. <|69|) ') can be included in. Maybe in a system of binary 
pulsars, the higher order terms in Ea. (|69|l are important. For static metric and in 1PN level Formula Q72|l agrees 
to the known formula in |27j . 

We proceed by evaluating the remaining terms in i|72[l . From Eq. I|60|l and Eq. (|61|l . we find 



Gb 
1 



We also have 

l-F(B) 



iGa 



w{B) 
w(A) 



w 2 {B) A Wl (B)v B 3w(B)v 



2c 2 



2c 4 



v\ w 2 (A) iwi{A)vi 



2c 4 
bw(A)v 



2c 2 



k B • vb 

c\k B \ 
4 Wl {B) dx i 
c 4 ~dx~ 



2c 4 



1_ ( k B ■ v B 

72 



2c 4 



8c 4 

M 

8c 4 



(73) 
(74) 



2w(B) kfl'Vg 1 fk B -v B 
c 3 \k B \ c 



k B • w B t Aw(B) ( k B ■ v B 
\k B \ 



\ks\ 

1_ ( k B • \ B 



\k 



B 



(75) 



B \ K B\ C- 

Substituting Eas. l|63|l . I|64|) . I|73l) . (|74|l and (|75|l into Ea. (|72|l . we finally have a general formula for the comparison 
of clock rates in the solar system, on the 2PN level of precision: 



Are 



w{A) - w(B 



1 



fk A -VA kg ■ v B 



1_ /" k B • v B 

72 



2c 2 

k A 'Va k B • v B 



1 



\k B \ 

w(B) -w{A) 



c\k A \ c\k B \ 

k B • v B , k A ■ v A 



c\k B \ 



c\kA\ 



\ c\kA\ 



2w{A) 



c\k B \ 



c\k B \ c\kA\ 
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\ (\Vw{A) - ^ _ 1 + 6w(B) ^ v 2 b 



4 - 2«i4 - v B ) +4[MB)v B - Wi(A)vh - Wi (B) t^-t — 777— p- 



-Kw)(w)( w(A)+w(i?) ) + (w) (w _ w , 

+ 5/;^ + 0(5). (76) 

Ea. (|76[) gives the frequency shift from A to B to the order 1/c 4 . If we really want to apply this equation 
in the case of the solar system, w(A), w(B), Wi(A) and Wi(B) have to be expanded by means of internal and 
external multiple moments |28) . As the calculation is rather long, we may publish it in a separate paper. We 
have checked that the leading terms of Ref. are included in the expansion of Eq. 176|) . 

To sum up our discussion we can say that Formula l|76ll 'contains' the Doppler effect, transverse Doppler 
effect (relativistic Doppler effect), gravitational redshift and their complete coupling effects to 2PN level in the 
solar system. In addition there is a term which is the integral of the rates of change of the scalar potential along 
the null geodetic line from source A to receiver B. This is probably the most interesting result in our paper. 
Hopefully this integral term and the coupling effects can be tested in the future with a deep space explorer and 
are confirmed. 



5 Conclusions 

1. We have synthesized all known effects for the comparison of clock rates in one formula Eq.J37|). The 
synthesized formula contains additional coupling terms and a new integral terms and thus gives essential new 
but untested information. Therefore, this synthesized formula is not an end in itself, but a starting point for 
the further test work. We hope that, this work could contribute to the further comparison of clock rates, such 
as ACES mission planned in near future. 

2. The general form is valid for any metric gravitational theory, not only for general relativity. If we 
substitute the parametrized 2PN metric into the formula, Ea. (|76|l could include parameters. 

3. In the case of cosmology, the general form can be used for any linearly perturbed metric, in particular it 
allows to include the general Sachs- Wolfe effects. 

4. In section IV, we have discussed the comparison of two clock rates both on the earth and a space station 
with 2PN precision. In fact, the calculation of the higher precision might be done in a similar way, if we know 
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the metric to higher order. 

5. Eq.J23 or Ea. 1(72(1 say that the clock rates depend on the trajectory of the transmitted signal and on the 
metric (especially the scalar potential) varying with time. While i(27() is a general formula, Eq. 1(721) is valid only 
in the solar system; however if replace Ea. 1)70(1 by the integral of Ea. 1(69(1 . a general 2PN formula results. 

6. The general form may be taken as the basis for a starting point to compare clock rates at any two different 
space time points. For example, the frequency shift caused by gravito-magnetic effect (or Lens-Thirring effect) 
can also be considered in our scheme. 

Acknowledgments 

This work was supported by the National Natural Science Foundation of China (Grant No. 10273008). In the 
workshop "Relativistic Astrophysics and High-precision Astrodynamics" led by Prof. Wei-tou Ni, we had several 
very fruitful discussions, and especially we would like to thank Prof. Tianyi Huang for his useful discussion. 
We also would like to thank the referees for their comments and pointing out and correcting our mistakes. 

References 



[i 

[2 
[3 

K 

[5 

[6 

[7: 

[8 
[9 

[io; 

[11 



Holzwarth R et al. 2000 Phys. Rev. Lett. 85 2264 
Udem Th et al. 2001 Phys. Rev. Lett. 86 4996 

Salomon C and Veillet C 1996 in Proceeding of the Symposium on Space Station Utilisation, ESA-SP 385 
p.295 

Spallicci A et al. 1997 Class. Quantum Grav. 14 2971 
Ni W-T 2002 Lnt. J. Mod. Phys. D 11 (7) 947 

Ni W-T, Shiomi S and Liao A-C 2004 Class. Quantum Grav. 21 S641 
Linct B and Teyssandier P 2002 Phys. Rev. D 66 024045 

Poncin-Lafitte C, Linet B and Teyssandier P 2004 Class. Quantum Grav. 21 4463 



Alba D and Lusanna L 2003 gr-qc/0311058 



Alba D and Lusanna L 2005 gr-qc/0501090 



LISA study team 2000 LISA: a Cornerstone Mission for the Observation of Gravitational Wave ESA 
System and Technology Study Report ESA-SCI 11 July 2000 



[12] Hannay J H 1985 J. Phys. A: Math. Gen. 18 221 



15 



[13] Spallicci A 2004 Nuovo Cim. B 119 1215 

[14] Spallicci A, Morbidelli A and Metris G 2005 Nonlinearity 18 45 

[15] Berry M V 1984 Proc. R. Soc. A 392 45 

[16] Ives H E and Stilwell G R 1938 J. Opt. Soc. Amer. 28 215 

[17] Saathoff G et al. 2003 Phys. Rev. Lett. 91 190403 

[18] Vessot R F C et al. 1980 Phys. Rev. Lett. 45 2081 

[19] Rces M and Sciama D W 1968 Nature 217 511 

[20] Blumenthal G R et al. 1992 Astrophys. J. 388 234 

[21] Tully R B et al. 1992, Astrophys. J. 388 9 

[22] Mcszaros A 1994 Astrophys. J. 423 19 

[23] Tuluie R, Laguna P and Anninos P 1996 Astrophys. J. 463 15 
[24] Martinez- Gonzalez E, Sanz J L and Silk J 1990 Astrophys. J. 355, L5-L9 
[25] Weinberg S 1972 Gravitation and Cosmology (John Wiley & Sons, Inc.) 
[26] Xu C and Wu X 2003 Chin. Phys. Lett. 20 195 

[27] Adlcr R, Bazin M and Schiffer M 1975 Introduction to General Relativity (McGraw-Hill Book Co.) 
[28] Damour T, Soffel M and Xu C 1991 Phys. Rev. D 43 3273 

[29] Schneider P, Ehlers J and Falco E 1992 Gravitational Lenses (Springer, New York) 
[30] Gong Y-X and Xu C 2004 Chinese Astronomy and Astrophysics 28 1 

[31] Will C M 1993 Theory and experiment in gravitational physics (Cambridge University Press) p. 167 

[32] Damour T, Soffel M and Xu C 1992 Phys. Rev. D 45 1017 

[33] Damour T, Soffel M and Xu C 1993 Phys. Rev. D 47 3124 

[34] Damour T, Soffel M and Xu C 1994 Phys. Rev. D 49 618 

[35] Bcrtschinger E 1995 Cosmological Dynamics (Elsevier Sci. Pub.) 

[36] Sachs R K and Wolfe A M 1967 Astrophys. J. 147 73 

[37] Bardeen J M 1980 Phys. Rev. D 22 416 

16 



[38] Ellis G R F and Bruni M 1989 Phys. Rev. D 40 1804 

[39] Anninos P http://relativity.livingreviews.org/open?pubNo=lrr-2001-2fcpage=articlesn9.html 
[40] Birkinshaw M and Gull S F 1983 Nature 302 315 

[41] Birkinshaw M 1989 in lecture Notes in Physics 330 Gravitational Lenses (eds Moran J M, Hewitt J N and 
Lo K L, Springer- Verlag, Berlin) p. 59 



17 



